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Abstract 

The formalism of causal dynamical triangulations (CDT) provides us with a 
non-perturbatively defined model of quantum gravity, where the sum over histo- 
ries includes only causal space-time histories. Path integrals of CDT and their 
continuum limits have been studied in two, three and four dimensions. Here 
we investigate a generalization of the two-dimensional CDT model, where the 
causality constraint is partially lifted by introducing weighted branching points, 
and demonstrate that the system can be solved analytically in the genus-zero 
sector. 
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Introduction 

The idea of CDT, by which we mean the definition of quantum gravity theory via 
causal dynamical triangulations, is two-fold: firstly, inspired by earlier ideas in 
the continuum theory [1,2], we insist, starting from space-times with a Lorentzian 
signature, that only causal histories contribute to the quantum gravitational path 
integral. Secondly, we assume the presence of a global time-foliation. 

The formalism of dynamical triangulations (DT) provides a simple regulariza- 
tion of the sum over geometries by providing a grid of pieccwise linear geometries 
constructed from elementary building blocks (these are rf-dimensional simplices 
of identical size and shape if we want to construct a d-dimensional geometry, 
see [3, 4] for reviews). The ultraviolet cut-off is given by the edge length of the 
building blocks. The causal variant CDT also uses DT as the regularization of 
the path integral. A detailed description of which causal geometries are included 
in the grid can be found in references [5, 6]. 

We emphasize that the use of triangulations is merely a technical regulariza- 
tion of the assumed underlying continuum theory, in the same way a lattice can 
be used for regularizing a quantum field theory. By no means do we presuppose 
that space-time is literally made out of little simplices. Some support for the 
existence of an underlying (non-perturbative) continuum quantum field theory 
in higher dimensions has been provided in [6]- [11] and seems to be in qualitative 
agreement with independent analyses carried out using the renormalization group 
[12, 13, 14, 15, 16]. 

While the CDT model is defined as a sum over causal space-time histories 
(each with an appropriate weight), one can ask whether this causality constraint 
can be lifted. One motivation for introducing it was that unrestricted summation 
over space-time histories^ leads to a dominance of highly singular configurations 
in dimensions d > 2, which prevents the existence of a physically meaningful 
continuum limit of the regularized lattice theory [17, 18, 19]. On the other hand, if 
one takes the point of view that a maximal number of possible fiuctuations should 
be included in the path integral (while still leading to a meaningful result), one 
may wonder whether it is possible to reintroduce (a subclass of) configurations 
into the sum over geometries which correspond to metric structures with causality 
violations. The question we would like to pose is whether this can be done in a 
controlled manner - using the additional time-slicing structure present in CDT 
- which avoids the problems encountered previously by DT, corresponding to an 
unrestricted inclusion of all such configurations. 

Because of the ready availability of analytic tools and the existence of analyt- 
ical solutions, we will in a first step analyze the situation in two dimensions. In 
this context, the issue has l)een addressed previously in a 2d toy model, focussing 

^even when keeping the space-time topology fixed 
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on the effects of including a class of minimal wormholes in CDT, which can be 
said to violate causality only mildly and are much less abundant than general 
wormholes [20, 21, 22]. In the present work, we will look at the genus- zero sector 
of a generahzed model of two-dimensional CDT, which in principle allows for the 
inclusion of arbitrary space-time topologies, as well as "outgrowths" , that is, the 
sprouting of baby universes, and associates with them a weight depending on the 
gravitational coupling constant.^ A key observation is that requiring the propa- 
gator of the model to reduce to that of standard CDT when the bare coupling is 
taken to zero uniquely fixes the scaling of this coupling, leading to a continuum 
limit where branching processes occur, but are scarce compared to the situation 
in DT. 

The following sections deal with explaining this scaling argument, and with 
analytically computing the genus-0 propagator (or loop-loop amplitude) and cor- 
responding disc amplitude. The computation of higher-genus amplitudes in this 
framework is currently under way, and may open new perspectives on the issue 
of the sum over topologies in theories of quantum gravity, which are only appar- 
ent in a formulation that has at least some memory of the Lorentzian structure 
of space-time built in. However, we have as yet no definite statements to make 
about the properties of general higher-genus amplitudes, the summability of the 
genus expansion or a generalization of the model to higher dimensions. 

Two-dimensional CDT and Euclidean quantum gravity 

Two-dimensional quantum gravity is a wonderful playground for "quantum ge- 
ometry" , understood as the statistical sum over geometries. The reason for this 
is that the action is trivial as long as we ignore topology changes (and even then 
it is almost trivial). One can therefore use entirely geometric reasoning to de- 
rive relations between or properties of "Crcen's functions"^. In this context it is 
convenient to study the proper-time "propagator" , namely, the amplitude of ge- 
ometries with two space-like boundaries separated by a proper time (or geodesic 
distance) t. Although the proper-time propagator is a special amplitude, it has 
the virtue that other amplitudes, like the disc or cylinder amplitudes, can be 
calculated from it [29, 30, 31, 32, 5]. When the path integral representation of 
this propagator is defined in the Lorentzian domain, using CDT, we can asso- 
ciate with each of the causal, piecewise linear Lorentzian space-time geometries 
a unique Euclidean geometry. After this rotation we perform the sum over the 

^The first analysis of a CDT model with local "decorations" was made in [23]. 

^An early example is the proof [24] that the string tension of bosonic string theory (regular- 
ized using DT [25, 26]) does not scale, thus providing a simple geometric understanding of the 
impossibility of defining bosonic string theory in target space dimensions larger than or equal 
to 2. Other applications in non-critical string theory can be found in [27, 28]. 
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Euclidean geometries thus obtained. The sum is now different from the usual 
Euclidean sum over geometries, since it extends only over a strict subset of all 
Euclidean configurations, leading to an alternative quantization of 2d quantum 
gravity (CDT). In the end, we can rotate back the propagator from Euclidean to 
Lorentzian proper time if needed. In the remainder of this article we will stay in 
the Euclidean regime, as defined above. 

For ease of presentation, we will in the following use a continuum notation. A 
derivation of the continuum expressions from the regularized (lattice) expressions 
can be found in [5]. We will assume that space-time has the topology x [0, 1]. 
After rotation to Euchdean signature, the action is 

S[g^.u]^X J d'e^det g^^iO + xj>dh + yj>dh, (1) 

where A is the cosmological constant, x and y are two so-called boundary cos- 
mological constants, g^^ is the metric of a geometry of the kind described above, 
and the line integrals refer to the lengths of the in- and out-boundaries induced 
by g^jii,. The propagator Gx{x,y;t) is defined by 

Gx{x,y;t)= J I^b^.] e"^!^-!, (2) 

where the functional integration is over all causal geometries [gfiu] such that the 
final boundary with boundary cosmological constant y is separated a geodesic 
distance t from the initial boundary with boundary cosmological constant x.^ 
Calculating the path integral (2) with the help of the CDT regularization and 
taking the continuum limit as the side-length a of the simplices goes to zero leads 
to the equation [5] 

-Gx{x, y; t) = -— [{x' - X)Gx{x, y; t)\ , (3) 

which is solved by^ 

Gx{x, y; t) = \ \ , 4 

x^ — X x{t, X) + y 

where x{t,x) denotes the solution of the characteristic equation for (3), namely, 

^ = -(x^-A), x{0,x)^x. (5) 



"^This class of geometries is difficult to define directly in a continuum, gauge-fixed formula- 
tion; what we have in mind here is the continuum limit of the corresponding discrete sum. 

^The asymmetry between x and y is due to the convention that the initial boundary contains 
a marked point. Symmetric expressions where neither or both boundaries have marked points 



can be found in [33]. 
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Figure 1: Graphical representation of relation 8: differentiating the disc amplitude 
Wx{x) (represented by the entire figure) with respect to the cosmological constant A 
corresponds to marking a point somewhere inside the disc. This point has a geodesic 
distance t from the initial loop. Associated with the point one can identify a connected 
curve of length I, all of whose points also have a geodesic distance t to the initial loop. 
This loop can now be thought of as the curve along which the lower part of the figure 
(corresponding to the loop-loop propagator Gx{x,l;t)) is glued to the cap, which itself 
is the disc amplitude Wx{l). 

Let /i denote the length of the initial and I2 the length of the final boundary. 
Rather than considering a situation where the boundary cosmological constant 
X is fixed, we will take li as fixed, and denote the corresponding propagator by 
Gx{li,y]t), with similar definitions for Gx{x,l2'it) and Gx{li,l2','t)- All of them 
are related by Laplace transformations, for instance. 



Eq. (7) is the simplest example of the use of quantum geometry. While the 
property (7) is evident in the context of CDT where no baby universes are allowed, 
it is also true in Euclidean quantum gravity (where there is no such constraint), if 
one defines the distance between the initial and final loop appropriately [29, 34]. 

Another, slightly more complicated example is illustrated graphically by Fig. 
1, which implies the functional relation 





where the Laplace-transformed propagator obeys the composition rule 
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It encodes the following: let Wx{l) denote the disc amplitude, i.e. the Hartle- 
Hawking amplitude with a fixed boundary length /, and Wx{x) the corresponding 
Laplace-transformed amplitude where x is a fixed boundary cosmological con- 
stant. Differentiation with respect to the cosmological constant A means marking 
a point in the bulk, as shown in the figure. Each configuration appearing in the 
path integral has a unique decomposition into a cylinder of proper-time extension 
t, (where the proper time is defined as the geodesic distance of the marked point 
to the boundary), and the disc amplitude itself, as summarized in eq. (8). 

Starting from a regularized theory with a cut-off a, it was shown in [5] that 
there are two natural solutions to eq. (8). In one of them, the regularized disc 
amplitude diverges with the cut-off a and the geodesic distance t scales canonically 
with the lattice spacing a according to 

Wreg — > a^Wxix), T] < 0, (9) 
a— >0 

Ueg t/a^, e — 1. (10) 

a— >0 

In the other, the scaling goes as 

Wreg ^ coust . + a'' W^A (a:) , ri = 3/2 (11) 

a— >0 

treg > t/a', £=1/2, (12) 

a— >0 

where the subscript "re^' denotes the regularized quantities in the discrete lattice 
formulation. The first scahng (9)-(10), with rj — —1, is encountered in CDT, while 
the second scaling (11)-(12) is realized in Euclidean gravity, i.e. Liouville gravity 
or gravity defined from matrix models. 

As demonstrated in [5], it is possible to treat both models simultaneously. 
Allowing for the creation of baby universes during the "evolution" in proper time 
t (by construction, a process forbidden in CDT) leads to a generalization of (3), 
namely, 

a'-^Gx,g{x, y; t) = [(«(^' - A) + 2^ a'^-'Wx,g{x)') Gx,g{x, y; t)] , (13) 

where we have introduced a new coupling constant g, associated with the cre- 
ation of baby universes, and also made the additional dependence explicit in the 
amplitudes. In [5] it was noted that for g = 1, that is, viewing this creation as a 
purely geometric process^, one obtains Euclidean quantum gravity. This happens 
because according to (9) and (11), we have either rj = —1, which is inconsistent 



^By this we mean that each distinct geometry (distinct in the sense of EucHdean geometry) 
appears with equal weight in the sum over two-dimensional geometries. 
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Figure 2: In all four graphs, the geodesic distance from the final to the initial loop is 
given by t. Differentiating with respect to t leads to eq. (15). Shaded parts of graphs 
represent the full, g'^-dependent propagator and disc amplitude, and non-shaded parts 
the CDT propagator. 

with (13), or we have from (11) that rj = 3/2 and thus e = 1/2, which is consis- 
tent with (11). On the other hand, setting g = 0, thereby forbidding the creation 
of baby universes, leads of course back to (3). 

In a non-trivial extension of previous work, we will now allow for the possibility 
that the coupling g becomes a non-constant function g = g{a) of the cut-off a. A 
geometric interpretation of this assignment will be given in the discussion section 
below. Since we are interested in a theory which smoothly recovers CDT in the 
limit as — s> 0, it is natural to assume that t] = —1, like in CDT. Consequently, 
the only way to obtain a non-trivial consistent equation is to assume that g scales 
to zero with the cut-off a according to 



where gs is a coupling constant of mass dimension three, which is kept constant 
when a 0. With this choice, eq. (13) is turned into 



The graphical representation of eq. (15) (or (13) for (7 7^ 0) is shown in Fig. 2. 
Differentiating the integral equation corresponding to this figure with respect to 
the time t one obtains (15). The disc amplitude W\^g^{x) is at this stage unknown. 

Note that one could in principle have considered an a priori more general 
branching process, where more than one baby universe is allowed to sprout at 
any given time step t. However, one observes from the scaling relation (14) that 
the corresponding extra terms in relation (13) would be suppressed by higher 
orders of a and therefore play no role in the continuum limit. 




(14) 



di 



d_ - 
dx . 



{x^-\)+2g, iy,,,,(x))G'A,,.(x,i/;t)]. (15) 
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Figure 3: Graphical illustration of eq. (17). Shaded parts represent the full disc 
amplitude, unshaded parts the CDT disc amplitude and the CDT propagator. 



In the next section we will show that quantum geometry, in the sense defined 
above, together with the requirement of recovering standard CDT in the limit as 
Qs — > 0, uniquely determines the disc amplitude and thus G\^g^{x,y;t). 



The disc amplitude 

The disc amplitude of CDT was calculated in [5, 35]. In [5] it was determined di- 
rectly by integrating G\{li, I2 = 0; t) over all times. This decomposition is unique, 
since by assumption t is a global time and no baby universes can be created. In 
[35] it was shown that it could also be obtained from Euclidean quantum gravity 
(matrix model results) by peeling off baby universes in a systematic way. By 
either method one finds 

Wxix) = (16) 
X + V A 

for the disc amplitude as function of the boundary cosmological constant x. In 
the present, generalized case we allow for baby universes, leading to a graphical 
representation of the decomposition of the disc amplitude as shown in Fig. 3. It 
translates into the equation 



00 00 



Wx,,A^)=w['lix)+g,J dt J dhdk ih + l2)Gflix,h + k;t)W,,,MWx,,Ah) 



(17) 

for the full propagator W\^g^{x), where we have introduced a superscript (0) to 
indicate the CDT amplitudes, that is, 

wilix)^W,,g^=oix) = W,ix), (18) 
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and similarly for quantities which were defined in eqs. (16) and (4) respec- 

tively. The integrations in (17) can be performed, yielding 

Solving for Wx,g^ (x) we find 



[X 



2 



A) + W^A,.,(^) 



where we have defined 

Wx,,^ (x) = ^ (a;2 - A)2 + [gsWlJ^/X) +x-y/x). (21) 

The sign of the square root is fixed by the requirement that Wx^g^{x) — > Wx{x) for 
gs — > 0, and Wx,gX^) is determined up to the value Wx^gX"^)- We will now show 
that this value is also determined by consistency requirements of the quantum 
geometry. If we insert the solution (20) into eq. (15) we obtain 

-Gx,gSx,y;t) = -—^WxMGx,gMy;t)\. (22) 

In analogy with (4) and (5), this is solved by 

Gxg {x,y;t) = ^'-^f , \ , (23) 

where x{t^x) is the solution of the characteristic equation for (22), 

-^--Wx,gAx), x{0,x)^x, (24) 

such that 

dy 



(25) 

m Wx,gM 

Physically, we require that t can take values from to oo, as opposed to just in a 
finite interval. Prom expression (25) for t this is only possible if the polynomial 
under the square root in the defining equation (20) has a double zero, which fixes 
the function Wx,g^{x) to 

Wx,g. (x) ^(x-a) y/(x + ay - 2gs/a, (26) 

where 

a = uVx, -u + ^ ^0. (27) 



In order to have a physically acceptable Wx^g^{x), one has to choose the solution 
to the third-order equation which is closest to 1. Quite remarkably, one can also 
derive (26) from (20) by demanding that the inverse Laplace transform Wx^g^{l) 



fall off exponentially for large In this region W\^g^{x) equals W^g^{x) plus a 
convergent power series in the dimensionless coupling constant Qs/X^^^. 

One can check the consistency of the quantum geometry by noting that using 
(23) in (8) the integration can be performed to yield 



9^ Wx,M 

which is indeed satisfied by the solution (20). 



(28) 



The loop-loop amplitude 

We mentioned above that the loop-loop propagator can be regarded as a building 
block for other, more conventional "observables" in 2d quantum gravity. One of 
the most beautiful illustrations of this and at the same time a non-trivial example 
of what we have called quantum geometry is the calculation in 2d Euclidean quan- 
tum gravity of the loop- loop amplitude from the loop-loop proper-time propagator 
[30] . The full loop- loop amplitude is obtained by summing over all Euchdean 2d 
geometries with two boundaries, without any particular restriction on the bound- 
aries' mutual position. This amplitude was first calculated using matrix model 
techniques (for cylinder topology) [36]. 

To appreciate the underlying construction, consider a given geometry of cylin- 
drical topology. Its two boundaries will be separated by a geodesic distance t, in 
the sense of minimal distance of any point on the final loop to the initial loop. It 
follows that we can consider the geometry as composed of a cylinder where the 
entire final loop (i.e. each of its points) has a distance t from the initial one and 
a "cap" related to the disc amphtude, as illustrated in Fig. 4(a). One can now 
obtain the loop-loop amplitude by integrating over all t and all gluings of the cap 
(we refer to [30] for details). An intriguing aspect of the construction is that the 
decomposition of a given geometry into cylinders and caps is not unique. One 
can choose another decomposition consisting of two cylinders of length ti and t2, 
with ti + t2 = t, joined by a cap, as illustrated in Fig. 4(b). As shown in [30], the 
end result is indeed independent of this decomposition. 

The whole construction can be repeated for our new, generalized CDT model, 
in this way defining a loop-loop amplitude. More precisely, although an exact 
equality of amplitudes corresponding to different decompositions like those de- 
picted in Fig. 4(a) and (b) is not immediately obvious at the level of the triangula- 
tions of the discretized theory''', the continuum ansatz (29) below is self-consistent, 

''because of the different arrangements of the proper-time slicings 
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(a) 



Figure 4: Two different ways of decomposing the loop-loop amplitude into proper-time 
propagators and a disc amplitude. Two points touch in the disc amplitude W, pinching 
the boundary to a figure-8, which combinatorially implies a substitution Wx^g^{l) 
lWx^g^{l) in the formulas. The time variables are related by ti + 12 = t. 



in the sense that it leads to a non-trivial symmetric expression for the amplitude 
with a well-defined Qs —>■ limit. The algebra is similar to that of [30]. 

We will denote the loop-loop amplitude by Gyg^{x, y), and its Laplace trans- 
form by Gx,g^{h,h), related in the same way as was discussed for the loop- loop 
propagator (c.f. eq. (6) and the discussion leading up to it). The integral equation 
corresponding to Fig. 4 is given by 

poo /"OO 

Gx,gXhM)= dt dlGx,gAlul;t)lWx,gSl + l2)- (29) 
Jo Jo 

Laplace-transforming eq. (29), the integrals can be performed using eqs. (23)-(26). 
After some non-trivial algebra one obtains 

where we are using the notation 



fix) = ^{x + aY-2gJa = Wx,gX^)/{x - a). (31) 
In the limit Qs ^ ^ one finds 

^'^^^ 2v^(x + V^)2(y + v^)2' ^ ' 
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a result which could of course also have been obtained directly from (29) using (4), 
(5) and (16). We note that the corresponding expression in the case of Euclidean 
2d quantum gravity is given by 

G'f\x,y) = , , , , , , , , , h(x) = y X + a/A, (33) 

which can be obtained from expressions similar to (23)- (26), only with W\^gXx) 
replaced by the Euclidean disc amplitude 

W^^\x) = {x- VX/2) h{x). (34) 

We observe a structural similarity between (30) and (34), with the function f{x) 
having the same relation to Wx^g^ (x) as h(x) has to wl^\x). The existence of well- 
defined, symmetric expressions for the unrestricted loop-loop amplitudes in our 
generalized CDT model (at genus 0) and thus in standard two-dimensional CDT, 
formulas (30) and (32), gives strong support to the claims that (i) the proper- 
time propagator does indeed encode the complete information on the quantum- 
gravitational system, and (ii) following the arguments given in [30] concerning the 
decomposition invariance of the loop-loop amplitude (c.f. Fig. 4), the continuum 
theory is diffeomorphism-invariant. 



Discussion 

The generalized CDT model of 2d quantum gravity we have defined in this paper 
is a perturbative deformation of the original model in the sense that it has a 
convergent power expansion of the form 

oo 

WxA^) = J2''n{x,X)[^y (35) 

n=0 

in the dimensionless coupling constant Qs/X^^^. This implies in particular that 
the average number (n) of "causality violations" in a two-dimensional universe 
described by this model is finite, a property already observed in previous 2d 
models with topology change [20, 21, 22]. The expectation value of the number 
n of branchings can be computed according to 

(n) - dW^A,,.(x) 

which is finite as long as we are in the range of convergence of Wx^g^x). As 
already mentioned, this coincides precisely with the range where the function 
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W\^g^{x) behaves in a physically acceptable way, namely, Wx^g^{l) goes to zero 
like exponentially in terms of the length / of the boundary loop. The same is true 
for the other functions considered, namely, Gx^gXhih^i) and Gx^gXh-ih)- 

The behaviour (36) should be contrasted with that in 2d Euclidean quantum 
gravity, and is reflected in the different scahng behaviours (10) and (12) for the 
time t. These scaling relations show that the effective continuum "time unit" in 
Euclidean quantum gravity is much longer than in CDT, giving rise to infinitely 
many causality violations for a typical space-time history which appears in the 
path integral when the cut-off a is taken to zero. This phenomenon was discovered 
in the seminal paper [29] . 

As we have already mentioned in the introduction, the calculations presented 
here should be seen as pertaining to the genus-0 sector of a generalized CDT 
model, which also includes a sum over space-time topologies. Although we have 
not given a precise definition of the higher-genus amplitudes in this paper, one 
would expect them to be finite order by order. If the handles are as scarce as are 
the baby universes in the gcnus-0 amplitudes, it might even be that the sum over 
all genera is uniquely defined. Whether or not this is so will clearly also depend 
on the combinatorics of allowed handle configurations. 

In the context of higher-genus amplitudes, it is natural to associate each han- 
dle with a "string coupling constant" , because one may think of it as a process 
where (one-dimensional) space splits and joins again, albeit as a function of an 
intrinsic proper time, rather than the time of any embedding space. An explicit 
calculation reveals that in the generalized CDT model this process is related with 
a coupling constant gl [37], which one may think of as two separate factors of gs, 
associated with the splitting and joining respectively. 

How does the disc amplitude fit into this picture? From a purely Euclidean 
point of view all graphs appearing in Fig. 3 have the fixed topology of a disc. 
However, from a Lorentzian point of view, which comes with a notion of time, it 
is clear that the branching of a baby universe is associated with a change of the 
spatial topology, a singular process in a Lorentzian space-time [38] . One way of 
keeping track of this in a Wick-rotated, Euclidean picture is as follows. Since each 
time a baby universe branches off it also has to end somewhere, we may think 
of marking the resulting "tip" with a puncture. From a gravitational viewpoint, 
each new puncture corresponds to a topology change and receives a weight 1/Gn, 
where Gjv is Newton's constant, because it will lead to a change by precisely this 
amount in the two-dimensional (Euclidean) Einstein-Hilbert action 



Identifying the dimensionless coupling constant in eq. (13) with g{a) = e 




(37) 
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one can introduce a renormalized gravitational coupling constant by 

TT— = tAt + - In Aa^. (38) 

This implies that the bare gravitational coupling constant Gisf{a) goes to zero 
like 1/1 In I when the cut-off vanishes, a — > 0, in such a way that the product 
ei/G'jv"/A^/^ is independent of the cut-off a. We can now identify 

as the genuine coupling parameter in which we expand. 

This renormalization of the gravitational (or string) coupling constant is rem- 
iniscent of the famous double-scaling limit in non-critical string theory^. In that 
case one also has Qs oc 6"^/*^^^"^ the only difference being that relation (38) is 
changed to 

^ ^ +^lnAa2, (40) 



GT G^ia) 4 

whence the partition function of non-critical string theory appears precisely as a 
function of the dimensionless coupling constant Qs/X^^^- 
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